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Summary 


Lyapunov stability theory is used to analyze the stability properties 
of equilibrium solutions to some linear and nonlinear partial differential 
equations. Stability conditions are established for systems which have, as 
a linear part, a self-adjoint differential operator or one which can be 
transformed into a self-adjoint differential operator. 

Although the techniques described are applied to systems occurring in 
hydrodynamics the methods can be used for similar problems in other fields 


as well. 



1. Introduction 


Frequently the analysis of systems describing physical processes gives 
rise to a stability problem of a system of partial differential equations. 
Often a stability analysis is carried out on an approximate system model 
having a finite number of degrees of freedom, usually obtained by a spatial 
discretization or a modal truncation method. The stability conditions so 
derived are sometimes not sufficient for stability except in the case of 
infinitesimally small perturbations. 

This approach is followed in one of the most recent contributions to 
the stability theory of hydrodynamical systems by Eckhaus . His theory 
for analyzing the stability properties of the solutions to nonlinear partial 
differential equations is based on asymptotic expansions with respect to 
suitably defined small parameters and series expansions in terms of the 
eigenfunctions. This method becomes cumbersome for more complex systems. 

Some of the systems studied by Eckhaus ^ are investigated here using 

r 2 31 

general Lyapunov stability theory. L ’ The context of the approach is 
the same as in ^ , that is, an investigation of the formal properties of 
certain mathematical relationships without a rigorous justification of these 
formal properties. It turns out, for the examples cited, that considerable 
improvement in stability conditions is obtained by the procedure recommended 
in this report, even to the extent of allowing nonlinearities. 

The stability of the equilibrium solution is defined in terms of the 
norm induced by the inner product of the Hilbert space on which the solutions 
of the system are defined. For a certain class of differential operators 
the stability conditions of the equilibrium solutions can be derived in a 
straightforward way. This class of differential operators can be considerably 
extended by a transformation of the differential operator which involves a 
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modification of the inner product. The general approach, even to systems 
containing nonlinearities, is illustrated by examples. 

2 . Statement of the Problem 

Many physical systems are formally described by partial differential 
equations of the form 

9 ii(t,jc) 

g— ^ — +Lu(t,x) = _0 (xefl; te [0, 00 )) (2.1) 

where u/t,x.) is an n-vector function and L is a matrix whose elements 
are linear or nonlinear differential operators specified on a bounded 
connected open subset of an m-dimensional Euclidean space, E™. The 
parameters of _L can be space dependent but not time dependent. In order 
to uniquely specify solutions of (2.1) a set of additional constraints or 
boundary conditons must be given, generally by a relation of the form 

Hu (t,x') =0 (x* e 9fi; t e[0,~)) (2.2) 

where H is a matrix whose elements are formally specified differential 
operators and 9fi is the boundary of Cl. Furthermore any solution will 
depend on some initial function _0^ (x.) belonging to the n-dimensional 
space of functions, 0 which we will assume is a Hilbert space with elements 
smooth enough to assure that solutions to (2.1) and (2.2) exist and belong 
to 0. 

A solution to (2.1) and (2.2) will be designated as ju(t, 2 c; _6 q) , 
that is, the solution starting at t=0 and with initial condition _0„ (x) e 0, 

i(°» J*q) = T be solution of particular interest is the equilibrium 

solution u (x) e 0, which is assumed to be u, = 0_. 
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In the following sections the stability of the equilibrium solution 
or the trivial solution, ju = _0, will be formulated and for some particular 
systems further investigated. 

To formulate the concept of stability we suppose that 0, as a Hilbert 

space, has an inner product <v,u> and the norm, ||.||> induced by the 

1/2 

inner product | J uj | = ( < ti,u i > ) - Then the distance at any time between 

and u^ ® is given by | |_u^ - | | . As used here a general inner 

product of two functions u, v e 0 is denoted by: 

CV,JJ> = / v T W(x) u df2 (2 . 3) 

T 

where ( ) denotes the transpose, W (jc) is a "weighting" matrix. Its elements 
can be chosen as continuous functions in x. such that 

W (x) = W T (x) (2. A) 

and 

6^ ja T u, > u T W(x) u. > $2 ^ “ > B-j > > ® (2.5) 

for all x. e ^ = ^ +3 ^‘ 

The conditions (2.4) and (2.5) assure that the norm induced by (2.3) and 

the norm induced by <v, _u>=/v T udQ are equivalent. 

Q 

As will be shown, an advantageous choice of W(x) depends on the form 
of the linear part of the operator L in (2.1). 

3. Stability of the Equilibrium Solution 

The concept of stability can be defined in many different ways. Here 
stability will refer to stability in the sense of Lyapunov, that is, a system 
is stable if for a sufficiently small perturbation from equilibrium, the 
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solution will remain close to equilibrium for all future time. The 
stability of the trivial solution ij = 0. can now more precisely be defined 
in terms of the norm given on 0. 

Definition . The trivial solution ia = C) of (2.1) is said to be stable in the 
sense of Lyapunov if for every real number e>0, there exists a real number 
6>0 such that I |^q(x) I I < <5 implies ||u,( t, x.5 ) | |< e for all t > 0. 

Definition . The trivial solution _u = 0, of (2.1) is said to be asymptotically 
stable if it is stable and in addition I |u.(t»x.j .0^) | |-H1 as t-**>. 

The stability properties of the trivial solution ja = 0_ of (2.1) can 
often be determined by consideration of the rate of change of the functional: 

V(u) = | |uj | 2 = < u.)U_> - / u J W.(x) .u d (3.1) 

a 

The time variation of V (u? along solutions is given by 

V(t) = V(u(t, x; ^(x)) t e [0,-) 

with V (0) = V ( (x) ) 

and its time derivative by 

= lim [ i (V(t+h) - V(t))] (3.2) 

dt h+0 h 

for ti (t,x; (x)) if the limit exists. 

Formally, from (3.2), it follows that since 

<ii(t + h) , _u(t + h) > - < u (t), u. (t)> = 

= <ia(t + h) + jj (t) , u (t + h) - ii(t)> , (3.3) 
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^ M t<u(t + h) + u( t ) , * (t + y . - ^ >-1 


(3. A) 


or 


dV(t) 

dt 


-2<«ct, x; ja (x) > 


yielding 


^ - 2 <u,L u> 


1 ) 


It is clear that if can be reduced to the form: 

dt 


(3.5) 


2HSML < 2a ||u(t,x; 0^) | | 2 = 2a V(u) 
dt ^ 


(3.6) 


where a is some real constant, it follows by integration that 


V ( _u) < V(^q) exp 2at 

or | |u (t, x; Iq)!! < IIOq^II exp “t (3.7) 

and the equilibrium solution u, = 0, will be asymptotically stable for a<0 
and stable for a < 0 for all t > 0. The functional V (u) with these 

r 3 ] 

properties is a Lyapunov Functional 1 . Thus, with the proper choice of 

inner product or norm, the square of the norm becomes the Lyapunov 

functional which establishes asymptotic stability. 

In the following section will be evaluated for a certain class 

dt 

of linear differential operators JL. The section is followed by examples in 
which these methods are extended to certain nonlinear systems. 

A . Self-Adjoint Linear Differential Operators 

Consider the system as given by (2.1) and (2.2) and let 1 be a linear 
operator. The extent to which one is able to reduce to the form 

^ In a completely functional analytic presentation, this equation would 

only be defined on the domain of L^, D(L) £ 0, but provided this domain is 
dense in 0, the conditions derived later are sufficient for stability 1 . 
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(3.6) depends intrinsically on the operator L. With respect to an inner 
product, L*, the adjoint operator of ^ is formally defined by 

<v, L u, > = <L.* v, u> (A . 1) 

for all ii in the domain of L. An important class of operators is formed by 
the self-adjoint operators. An operator is formally self-adjoint if L* = L. 

If one restricts oneself to the inner product 

<v,_u> = / v T ii dfi (A. 2) 

n 

it is clear that the corresponding class of self-adjoint operators is very 
limited. The introduction of the matrix, W(x) , often permits to be made 
self-adjoint with respect to the generalized inner product (2.3). This is shown 
in the following example. 

Example . Let the differential operator be given by 

L u = a (x) + b (x) + c(x) u, (A. 3) 

0 £ x < 1, a(x) > 6 > 0, b(x) £ 0 and boundary conditions u(0) = u(l) = 0. 

Then with the inner product (A. 2): 

<v, L u> = v L u dx = {a(x)v ^ + b(x)v “ + c(x) vu}dx. 

0 0 

After integration by parts and substitution of the boundary conditons, 
assuming v(0) “ v(l) = 0 one gets 

1 j\2 g 

<v, L u> = / (u -^ 2 (a(x)v)-u (b(x)v) + c(x)vu}dx 

= <L*v, u> 

where L* is given by 
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(a(x)v) - — • (b(x)v) + c(x)v 



0 < x < 1 and boundary conditions v(0) ■ v(l) = 0. Clearly L* ^ L and L 
is not self-adjoint. 

But L is equivalent to L e as given by 


L u = 
e 


1 

w(x) 


8(p(x) fj> 
3x 


+ q(x) u. 


(4.4) 


0 < x < 1, where q(x) = c(x) 

p(x) = exp / dx w(x) = (a(x)) 1 exp / x dx 

and boundary conditons u(0) = u(l) = 0. This can be verified directly. 
Application of the general inner product (2.3) gives 

<v,L u> = v L u w dx = / (v (p(x) ~) + vq(x)uw)dx 
e 0 e 0 dx 

Integration by parts of the expression in the integrand and substitution of 
the boundary conditions, with the assumption v(0) = v(l) = 0 gives 


<v, L o u> = /■*■ (u ~ (p(x) ~) + u q (x) v w} dx = <L*^v, u> = 


3x 


= <L v, u >. 
e 

Thus will be self-adjoint with respect to the inner product 


<v, u>-/^v u w dx. 
0 


(4.5) 
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The significance of the choice of inner product becomes apparent on 
evaluation of given by (3.5) where V(u) is given by (3.1), for 

self-adjoint operators. If L is a self-adjoint operator with a lower semi- 
bounded spectrum and is the smallest eigenvalue of I., then 


^min 


, Ji> < <U., Ij u> 


( 4 . 6 ) 


To show how this property can be used, consider the linear system 


3ii 

r~ + L u = 0. (4.7) 


and let L be a differential operator with possibly space dependent coefficients. 
Let a set of boundary conditions be specified and let I. be self-adjoint with 
respect to the general inner product 

<v_, u> = / W(x) u dfi (4.8) 

n 

From (3.5) and the fact that _L is self-adjoint with a lower semi-bounded spectrum 

and X . is the smallest eigenvalue of L, it follows that 
min — 


dV(u) 

dt 


< -?.A . <u, u> 

= min — 


(4.9) 


Hence if X . > 0 we have (3.6) and thus asymptotic stability, 

min 

A similar approach can be followed when L is time varying by introducing 

W(x, t) . The elements of W^(x, t) must be continuous in both jc and t and 

continuously differentiable in t and such that (2.4) and (2.5) are satisfied. 

The derivative dV(ji) can be evaluated as in Section 3. This yields an 
dt 

additional term in (4.9) which can be estimated to yield a sufficient condition 
for stability. 
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The use of this general inner product has the advantage of relating 

the stability properties directly to the eigenvalues of the differential 

operator for a large class of operators. As will be shown in the examples 

in the next section in order to reduce ^ to the form a <u^ one 

dt 

can apply well-known integral inequalities rather than calculating the 

eigenvalues. 7116 use t,le general inner product representation facilitates 

the application of these inequalities and improves on the resulting stability 

[4] 

condition. Of course these inequalities , although standard, are to some 
extent based on estimates of the eigenvalues. 

When applied to nonlinear differential operators the advantages become 
even more apparent since this permits an estimate for the set of initial 
functions in 0 for which the trivial solution is asymptotically stable. In 
the next section this will be illustrated by applying the above methods to Burgers' 
model of turbulence and some of its modifications as studied by Eckhaus^^. 

5. Applications . 

Example 1. Burgers' model to describe turbulence as studied by Eckhaus ^ is 
given by 



1 

R 


9 2 u^ 9u^ 

9 x^ 9x 


U 1 U 2 


0 




dx = 


n 


(5.1) 


0 < x < 1 and boundary conditions u^(0) = u^(l) = 0. To illustrate the 

foregoing we will make a slight generalization of the above problem and 
in the sequel establish the results of Eckhaus. 
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(a) Consider first the linear system 


3u_ 

at" + iit = 2 . 





. 1 3 2 , . 


U 1 


- 1 - R 3TF- J 

where u = 


; L = 



U 2. 


° ’ R . 


with R = R(x) > 0 for xe [0,1] and boundary conditions u^ (0) = u^(l) = 0. 
This assumption complicates the problem, because with constant R the operator 
Ij is self-adjoint with respect to the inner product 

1 

<v, u> = / 0 (v 1 u 1 + v 2 u 2 )dx (5.3) 


and the condition for asymptotic stability of the trivial solution ii = Q_ 
can be derived from the smallest eigenvalue of _L which is easilj determined. 

However, since R(x) is assumed to be space dependent, some sort of 
variational technique would be required to obtain . This can be avoided as 

shown below by making use of an integral inequality. 

For the system as given by (5.2) choose as V ( ti ) functional: 



V(u) = | |u| 

rT-i 

| 2 = <ii, u> = f q u W(x) u_ dx 

(5.4) 

where 

W(x) 

_[r(x) . °1 

-[ o , xj • 

(5.6) 


V ( ) is now the required Lyapunov Functional ^ if — ^7^ 2- <0. is 
self-adjoint for the inner product (5.4), however determination of the 
eigenvalues is not as immediate as in the case of constant R. The derivative of 



dV(u) 


(5.7) 


— 7 — = -2 <u, L u> 
at — 

Substituting and W(x) gives: 


dV(u) 

dt~ 


8 z u, 


= -2 J 0 (-R(x)u l — u 2 }dx. 


Integration by parts of the second term in the integrand of (5.8) and 
tution of the boundary conditions gives 


dV(u) 


1 „ 3u „ 1 „ 

= -2 J q (-R(x) u L + (— ) + ^ up dx 


r 4 1 

Now the following integral inequality' holds for u 


1 3u 


/ ( ) 2 d x > n 2 / uif dx. 

J 0 = 0 1 


1 2 


Thus 


dV (u) 


dt <- 2 1 0 {(t 2 - r W )u i + rw “2 


u 2 }dx 


dV(u) 1 1 

at - i ~^ n [0>1] <7r2 “ R(x)) * roo 5 


The equilibrium solution ti = () will be asymptotically stable if 


min , 1 

xs [0,1] l R(x) 


(tt 2 - R(x)),^y ) > 0 


Since R(x) > 0 for x e [0,1] condition (5.13) reduces to 

it 2 > max R(x) • 

xe [ 0 ,l] 


(5.8) 
substi- 

(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 

(5.14) 
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With constant R>0, for u_ = 0_ to be asymptotically stable, it follows immediately 
from (5.14) that a sufficient condition is that 

0 < R < it 2 (5.15) 

This condition is identical to that obtained by requiring that X ° f L. be 
positive, since it is shown in (1] that. 

IT 2 

\ ln = ~ (n + 1) - 1 (n=0,l,2, . . .) 

» 20 -i (5.16) 

3u 2 

(b) Next consider the nonlinear case and suppress (as in [1]) the r~ 

O t 

term in (5.1). This gives: 


9u 1 

3t “ U R 



0 


(5.17) 


0 < x < 1 and boundary conditions u(0) = u(l) = 0. For simplicity, assume 
in addition that R is a positive constant here, ^aking as V(u) functional 


V (u) 



dx 


(5.18) 


its time derivative becomes: 


dV (u) „ r 1 r 2 1 , 3 2 u v , 3u 2 . „ , r 1 2, , 2, , 

— = -2 J { _u - - u(-^ 2 -) + u — + R ( /_ u dx) u }dx. 


R 3x^ 


3x 


(5.19) 


Integration by parts and substitution of the boundary conditions gives: 

dV (u) „ r 1 r 2 , 1 , 3u ,2 . „ rf l 2 . s 2, , , r 

— = -2 J l-u + — (— ^ — ) + R (J n u dx) a }dx (5.20) 


r 2 

Applying the integral inequality (5.10), the fact that u dx > 0 and 

R > 0 gives: 


^ s £ < r - « » 2 


dx. 


(5.21) 
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Thus the modified nonlinear system also has an asymptotically stable 
equilibrium solution u =* 0 for 


0 < R < it 2 


(5.22) 


This verifies the previously obtained result in [1] suppressing the 


8u 2 

TE“ 


term. 


(c) Next it will be shown that a stability analysis of the nonlinear system (5.1) 
can be made with the techniques described above without making the modification 
under (b) , a result not obtained in [1]. Consider the system as given by (5.1) 
and let R be constant. Since R is constant take as V(u) functional: 


1 , 2 


V(u) = < u_ , u_ > = (u^ + u 


;) dx 


(5.23) 


The time derivative of V(u) becomes after substituting (5.1): 


dV(u) . r 1 r 2 1 3 2 ui 2 .12. r 1 2 , . 3uf, , /r 

dT = ' 2 to [ " u i “ r u i ~a^T ~ u l u 2 + R U 2 + u 2 J 0 u i dx + U 1 lt ] dx (5 ' 24) 


Integration by parts and substitution of the boundary conditions and inequality 
(5.10) gives: 


= " 2 [( r” -1) u i + ¥ u 2 " u 2 (u r f 0 u ? dx)]dx * (5 - 25) 


The inequality (5.25) can also by written as 
dV (u) 


dt i " 2 /q [( r~ 1 )u l + R U 2 ]dx 


and certainly: 


1- 


f o u 2 (u r /J u i dx)dx 


Jo [< r- 1 )u i + r u 2 ] dx 


(5.26) 


dV(u) 

dt” 


J-2 / l 


1- 


l Q l u 2 (u l + jp u i dx) l dx 


(1 r f TT Z - v 2 1 2, , 

/ - 1 > u l + R U 2 ]dx 


(5.27) 
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Repeated application of the Buniakovsky — Schwartz inequality and the fact 
that u^dx < | |uj j 2 and /^u^dx < jjujj 2 yields from (5.27): 


dV (u) 


dt 


1 w 2 2 1 2 2 1 llil I 

-2 / [( ~— -l)uf + £ u^]dx 1- -2 i (5.28) 

0 R 1 R 2 «in(f--l # 1) ||u||2 


or 


dV (ji ) 
“ 


-2 / 0 [<ir-i) 


u l + l u 2 ]dX fl- I lull]- 

l rain( T J 


(5.29) 


A sufficient condition for the trivial solution u=0 of (5.1) to be 
asymptotically stable is thus that 

0 < R < TT 2 


(5.30) 


and the disturbances be bounded in norm by 

| |u| I <\ min (|- - 1 , i) . (5.31) 

Notice that (5.30 - 5.31) requires the linear approximation of the system (5.1) 
to be asymptotically stable, which should be expected. 

The conditions (5.30 - 5.31) are derived without any prior knowledge about 
the solutions, an important advantage of the Lyapunov approach. However an 
examination of (5.1) shows that ^ is independent of x. Consequently can be 
taken outside the integral sign in (5.26). Upon integration there follows 


dV (u) 


dt 


< -2 fl [( Jr -1) u? + i uh dx (5.26a) 

= '0 R 1 R 2 


A sufficient condition for the asymptotic stability of the trivial solution 
ii *= 0_ of (5.1) is thus 0 < R < tt 2 , the same as the linear approximation case. 
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This example shows that the Lyapunov stability theory not only enables us 
to verify the results of^ but to extend these results considerably. 

Example 2 . Next we consider a second example from ^ 


3u 

3t 


(x 


2. 

2 3u 

1 

8u 2 2 

l T) 

f 1 u 2 dxi 

f 0 U 

,/R* X 3x “ 

R 

™ "T K 

ax 

J 0 u d *J 

► 


u = 0 


(5.32) 


0 < x < 1, boundary conditions u(0) = u(l) = 0 and R a positive constant. 
The linearized system is given by: 


3u 

at 


(x 



2 3u 1 a 2 u 


(5.33) 


0 < x f 1 and boundary conditions u(0) = u(l) = 0. The linear differential 
operator L of (5.33) is not self-adjoint for the inner product 

1 

<v,u> = vu dx. (5,34) 


Hence on the basis of the preceding, we might expect that the stability 
condition derived with this inner product would not be the best possible one. 
To illustrate this take as the functional V(u): 

1 ? 

V(u) = J Q u dx. 


The time derivative follows after integration by parts as 


dV (u) 
dt 


„ r f l ,8u.2 ,1.2. 2, 

= - 2 Jo % ( a^ } - ( x > u }dx * 

Applying the inequality (5.10) gives 


(5.36) 


dV (u) 
dt 


- fl , li 1 2, 2 , 

“ 2 J 0 ( R "ft- X )U dX * 


(3.37) 
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And a sufficient condition for asymptotic stability of the solution u=0 is: 


or 



(x 2 }> 0 


(5.38) 


0 < R < j (1 + 2n 2 -\jl + 


(5.39) 


One can improve greatly on condition (5.39) as an evaluation of the 
eigenvalues of L suggests (see [1]) by observing that L is equivalent to the 
operator as given by: 


L u 
e 


w(x) 3x 


f^(p(*> |^J+ qU)u 


with 


p = exp 




w = 



2 

x 


2 2 
q = x 

and boundary conditions u(0) =■ u(l) = 0. 
The functional V(u) will be taken as 


V(u) 


< u,u 



2 

u 


w(x)dx. 


(5.40) 


(5.41) 


After integration by parts and substitution of the boundary conditions 
its time derivatives becomes: 

dV(u) . „ r r V/k x ,Su N 2 „ Jr x , 2 , 2, 2,j , c 

d r = - 2 Jo {e ( 3^ - Re (x + ^ )u }dx - (5 * 42) 

i^x 2 

Here the inequality (5.10) can again be applied, now however to e u, 

rather than to u. Substituting the result into (5.42) gives: 



dV (u) 
dt 


(5.43) 


± “ 2 /o ( ^ W(x) dx * 

The condition for asymptotic stability of the equilibrium solution u=0 of 
(5.33) follows from (5.43) as 

0 < R < tt* 4 . (5.44) 


This condition is identical to that found by evaluating the eigenvalues of L 

as should be expected. Further, comparison of (5.39) and (5.44) shows the 

intrinsic dependence of the stability condition on the V(u) functional chosen. 
The modification of the functional V(u) based on a transformation of the 

differential operator L in this case results in a significant improvement of 

the bounds of the system parameters to assure stability. However, as indicated 

in Example 1, the results are even more important when dealing with nonlinear 

systems. To show this consider the nonlinear system (5.32). Take as V(u) 

functional (5.41), thus 

V(u) = f u w(x) dx. (5.45) 

J 0 


The time derivative of V(u) becomes: 


<*V(u) = _ 2 f 1 {e x 2 ( ^u. 2 
dt 3x 





u 





+ 


t J u dx]u }dx. (5.46) 


Integration by parts, substitution of the boundary conditions and the inequality 

, 1 2 

(5.10) together with the fact } u dx > 0 and R > 0 give 

0 


dV(u) 

dt 


< -2 J 0 t ~ ~ J x u]u 2 w(x) dx. 


(5.47) 
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will certainly be negative definite for 


Since 0 < x < 


i dV<u) 

» dt 


R 



> 0 


(5.48) 


for all x e [0,1] . 


Hence the equilibrium solution u=0 of (5.32) will be asymptotically stable 
for all disturbances bounded by 


max 

xe [0,1] 



(5.49) 


This result is again similar to that obtained by Eckhaus [1], however 
the above procedure enables one to obtain it in a straightforward way without 
making many complicated calculations as is the case when using asymptotic 
expansions . 


6 . Conclusions 

Some recent results in the stability theory of partial different ial equa- 
tions have been obtained by means of Lyapunov stability theory. Although the 
application of this theory might not be as unified as the approach outlined 
by Eckhaus [1] it is felt that with some sophistication many results can be 
obtained in a less cumbersome way. It is also gratifying that many of the 
results of Eckhaus are verified using Lyapunov stability theory. Among the 
useful tools that can be applied we have demonstrated a transformation of the 
differential operators and the use of integral inequalities. 

The feasibility of Lyapunov stability theory in the analysis of solutions 
to partial differential equations is clearly established by the ease with 
which sufficient conditions for stability are determined for the linear and 
nonlinear systems discussed in the examples. 
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